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1 Introduction 



Two-dimensional Hydrogen atoms in magnetic fields have been a subject of active 
research during the last years. A large body of articles has been published on this 
problem in the framework of non-relativistic quantum mechanics. This problem is 
of practical interest because of the technological advances in nanofabrication tech- 
nology that have made possible to create low dimensional structures like quantum 
wells, quantum wires and quantum dots |1], |2|, [J. The two-dimensional Hamiltonian 
describing the Coulomb interaction — — , between a conduction electron and donor 
impurity center when a constant magnetic B field is applied perpendicular to the 
plane of motion, can be written in atomic units, h = M = e = lin the CGS system, 
as follows 

11 Z 
H<p = -(-iV + -Bxf f ) 2 <p--<p = id t <p = E<p (1) 

Since we are dealing with a two-dimensional problem, we choose to work in polar co- 
ordinates (r, The angular operator operator — id& commutes with the Hamiltonian 
([!]), consequently we can introduce the following ansatz for the eigenfunction 

= expfomfl) u{r) ^ 
v 2tt 



Substituting (0) into (P, we readily obtain that the radial function u{r) satisfies the 
second order differential equation 



Id 2 1, 2 1,1 LU 2 L r 2 Z 



u{r) = 0, (3) 



where u>l = B /2c is the Larmor frequency, E is the energy, and m the eigenvalue of 
the angular momentum. Eq. (|3|) cannot be solved in closed form in terms of special 
functions 0]. There are analytic expressions for the energy for particular values of 
ojl and m[[5], || 0|. The computation of the energy eigenvalues in (|3]) has been carried 
out using different techniques |8|, ||, [10], [□], [L2|, pUf . A less studied problem is that 



of a relativistic 2D hydrogen atom in a magnetic field. Perhaps relativistic effects 



1 



are not very important in semi-conducting devices but nevertheless they cannot be 
neglected when the interacting potentials are strong [ 14] , |T5| . Recently, the importance 



of considering relativistic effects has been pointed out when one computes the energy 
levels of semiconductors in high magnetic fields |16|. The effective mass method is 



still valid until H « 10 5 Oe. In the case of two band approximation, the dispersion 
law has the form of a Klein-Gordon energy spectrum]!?]]. Also, relativistic invariance 
imposes some supplementary restrictions on the allowed quantum energy levels. 

In this article we investigate the relativistic corrections to the energy spectrum of a 
two-dimensional hydrogen atom in an homogeneous transverse magnetic field. Using 
a mixed-basis variational approach |18| , |l9 |. In Sec 2, we compute the relativistic 



energy spectrum of a 2D relativistic Klein-Gordon hydrogen atom. In Sec 3, we 
discuss the application of the 1/N expansion to our problem. In Sec 4, we compare 
the energy spectrum of the relativistic 2D hydrogen atom with that obtained in the 
nonrelativistic limit. Finally, we present the concluding remarks in Sec. 5. 

2 Relativistic Hydrogen Atom 

Since we are interested in discussing the relativistic corrections to the energy levels 
of the 2D Hydrogen atom, we proceed to solve the 2D Klein-Gordon equation. The 
results obtained after solving the Klein-Gordon equation apply to an electron without 
spin. The advantage of this approach |15| can be easily understood if we recall that 



the Schodinger equation does not take into account the spin of the electron ([]]) and 
then we can directly compare the relativistic and nonrelativistic energy spectra. 

The covariant generalization of the Klein-Gordon equation in the presence of elec- 
tromagnetic interactions takes the form [ |TB| , 



(g aP (V a - -A a ){V p - -A p ) - c 2 ) * = o, (4) 
where the contravariant metric tensor, and V Q is the covariant derivative. The 



2 



metric tensor g a p written in polar coordinates (t, r, $) takes the form: 

g a p = diag(-l, l,r 2 ), (5) 

and the vector potential A a associated with a 2D Coulomb potential and a constant 
magnetic field interaction is 

7 Br 2 

^ = (—,0,-—). (6) 
r 2 

From the above expression (||) for the vector potential A a it is straightforward to 
verify that the electric and magnetic fields satisfy the invariant relations 



F aP F aP = 2(B 2 - E 2 ) = 2(B 2 - (7) 

*F a/3 F a P = ^E-B = } (8) 

where F a ^ is the (2+1) electromagnetic field strength tensor. 

Expressions (|7|) and (|8]) tell us that in fact, A a is associated with a 2D Coulomb 
atom in a constant magnetic field perpendicular to the plane of the particle motion. 
The corresponding E and B can be written in polar coordinates as follows: 

E = -—e r , B = Be z . (9) 

Since the vector potential components do not depend on time or the angular variable 
we have that the wave function solution of the Klein-Gordon equation can 
be written as 



¥(r, tf, t) = expUmtf - Et), (10) 



where the function u{r) satisfies the second order differential equation 



^ + A-m- + S _rng_ c2 + g_W + 2gg, 
dr 2 \ r 2 c c 2 4 c 2 c 2 r 
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Eq. (|ll]) has the same form as Eq. (Q), therefore no exact solutions of Eq. (|TTD can 
be obtained in terms of special functions. 

In the present article we analyze the problem of computing the energy levels of 
the 2D relativistic Coulomb atom using a mixed-basis variational approach. In order 
to apply the variational method to our problem [EOj , we look for a trial wave function. 



Since Eq. (|TT|) reduces to the relativistic Hydrogen atom equation when uo L = 0, we 
can consider as a basis, for ojl « 1, the Hydrogen wave functions %. The solution 
of eq. (^j when ojl = is 



»//(/■) = D WII , -i^-%r { J^ +X/2) L{n p ,<l\ln^ — .2/ V - ~) (12) 
where D m ^ n is a normalization constant, L(a,b,x) are the Laguerre polynomials 



Z 2 

2 



2 • 



and E from 22 



is 



E = c 2 



1 + 



c 2 (n 



P 2 



+ \ m 



-,2 _ 7^\ 2 



-1/2 



(13) 



It is worth mentioning that the relation (13) makes sense only when 



m > 0, 



(14) 



a condition that forbids the existence of the s energy levels (m = 0), this is in fact 
a peculiarity of the relativistic Klein-Gordon solution, which is not present in the 
standard Schrodinger framework. 

Conversely, for large values of ul, a good trial basis is that of the spherical oscil- 
lator. In this case, the solution of eq. (|TT|) has the form 



5 (r) = C min e-^ 2 /2 p (|m| + i) L( H j(jLf ?) 



(15) 



and, in the high-field limit, the energy spectrum of a relativistic spinless particle in a 
constant magnetic field satisfies the relation 
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E 2 

— - c 2 = 2u L {2n + m + \m\ + 1) . (16) 

<r 

Among the advantages of considering a relativistic spinless electron is that we 
can easily compute the energy levels with the help of the mixed base variational 
approach. Nevertheless, we can easily see when the Klein Gordon equation gives a 
reasonably good value for the energy spectrum as compared to that obtained via the 
Dirac equation. 

The 2+1 Dirac equation E7I, EBJ in the presence of an external electromagnetic 



7^ + ^)+c)*=0. (17) 



field Afj, reads 

( V/ 

v dx^ c 

Since we are working in a two dimensional space, we can work in the following repre- 
sentation of the gamma matrices 

7 ° = za 3 , y = (T i ; ^ = a \ (18) 

Then, the Dirac spinor has only two components. Since the Dirac equation ( |17D 
expressed in the diagonal tetrad gauge commutes with the operators and — i-j^, 
the spinor \l/ can be written as 

tf(t,r) = -^exp(-iEt + iW)ip,(r) (19) 
v 2n 



^\ Mr) \ (20) 



ip 2 (r) 

where ip\ and ip2 satisfy the system of equations 

,E Z, ,d I Bp, 

c p op p 2c 

/ E Z s , . d I Bp. , 

( + c _ _ ^ + + ^ ^ = o. (22) 

c p op p 2c 

Substituting fl2l|) into (|2|), we obtain the following second order differential equation 



d 2 ^ 2 /,£ Z N9 9 1(1 + 1) lB 2 p 2 „ l,B 

+ (- + -f - c 2 - - -— f- -(/ + -)- + £> U 2 = (23) 



dp 2 \ c p p 2 A c 2 2 c 
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where D is given by 



„ Z ,E Z, u d I Bp.. 

£> = — (- + c+-) 1 ^- + - + ^V'2 24 
p 2 c p op p 2c 



It is worth mentioning that angular parameter I takes half integer values |27j and 
therefore it can be related to m as follows 

l = m--. (25) 



Eq. (p3[) reduces to the radial Klein-Gordon equation (|TTJ) when .D vanishes. 
Looking at (|24| ) and (|21"| ) we see that L> is very small for large values of p |[28|| . The 
mean square radius of the Dirac and Klein-Gordon electron is p 2 ~ 2(n + l/2)c/B 
where n labels the energy levels. Taking into account that large values of the radial 
variable imply that p > 1/c, we find that D is negligible for magnetic fields satisfying 
the inequality B < 2c 3 , which is the critical value for strong magnetic fields |29| . 



Using Eq. (0), and keeping only leading terms of D, the motion of a relativistic 

electron for small values of B and p is described by the equation 

(Pu> 2 ( ,E Z,~ 9 j -I 2 lB 2 p 2 R\ 

+ [(- + -) - c 2 + i - 5 - " T— H -(I- !)- ^2 = 0. (26) 



tip 2 \ c p P 2 4 c 2 c 

Notice that Eq ( PB[ ) gives a description of a relativistic electron in a weak magnetic 
field B for small values of p. From Eq. (p3|) and (|25|, we see that, in oposition to the 



Klein-Gordon case, the 2+1 relativistic Dirac electron has s states. 

If we attempt to apply the variational method using the hydrogen atom basis, 
we will obtain good agreement with accurate results for small values of cjl, but this 
approach fails for large u>l, even if we consider a many term basis. An analogous 
situation occurs when we use the oscillator basis, in which case we obtain a good 



agreement for large ujl, but the convergence is very slow for small values of uj^. | 



In order to solve this problem, we propose a trial function |T8], [19fl , for any quantum 
level n, a linear combination of the form 



N 

U 



^2(c iH UiH + CiQUiOsc) (27) 
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where N > i > n; itm and UiOsc are the corresponding hydrogen and oscillator 
wave functions associated with the quantum level % Cio and cm are constants to be 
calculated. It is worth mentioning that our basis is not orthogonal under the inner 
product (ui | Uj) = J °° UiUjdp. Substituting ( p7|) into flllD , and performing variation 
on the basis coefficients Cj, we readily obtain the following matrix equation: 




where the indices i and j running from 1 to N correspond to the Hydrogen and 
oscillator bases respectively. The algebraic equation QijCj = gives nontrivial 
values of Cj provided that the matrix Qij be singular. The energy eigenvalue for a 
given quantum level n is the lowest value of E, solution of the equation det(Qy) = 
0. The mixed-basis variational method gives reasonably good values for the energy 
eigenvalues even for the simple selection of a two term basis as in Eq. (|27|) . In this 
particular case, we have that the trial function, for any quantum level, is a linear 
combination of the form 

Ui = c iH u iH + c i0 u i0 (30) 

Better results should be expected for a basis with more terms. For a three term basis, 
we have two possible trial functions 

Ui = CiUuh + C 2 U 2 iH + C 3 Uuo, (7712X21) (31) 

and 

Ui = CiUuh + c 2 uii + c 3 u 2i o, (mixl2) (32) 



7 



where the three terms in (|31]) and (|32|) have the same angular dependences of the 
eigenfunction to be approximated. In this scheme U2« corresponds to a wavefunction 
associated with a higher quantum number to that we are going to approximate. 

For comparison, the numerical computations of the relativistic energy spectra are 
carried out with the help of the Schwartz method []24|] , which is a generalization of the 
mesh point technique for numerical approximation of functions. This method gives 
highly accurate results given a thoughtful choice of the reference function, and its 
efficiency has been shown computing the energy spectrum of the 2D Hydrogen atom 
0. 



3 1/N Approach 

The shifted 1/N expansion is a perturbative technique that has permitted us to solve 
the N-dimensional stationary Schrodinger equation with a wide class of radial poten- 
tials. The shifted 1/N method has also been developed to compute energy eigenvalues 
of relativistic spin and spin | particles in the presence of spherically symmetric vec- 
tor and scalar potentials. Here we proceed to compute the energy eigenvalues of 



our problem using the 1/N expansion for N = 2. Since Eq. (11) contains a mag- 



netic field contribution, some minor modifications should be made to the recipe of 



Ref [25[] where the authors develop the shifted 1/N technique to deal with the 3D 



Klein-Gordon equation in a spherically symmetric potential. 



Following the scheme developed by Imbo and Pagnamenta | ID| , we have that the 
radial Klein-Gordon equation in the presence of a constant magnetic field @ takes 
the form 

£ (k + a-l)(k + a-?>) , _ D , „ , 1 2d2 



, 2 + j\ + mB + M 2 + -r l B l - [E + -] 2 U nr (r) = 

ar l 4r 2 4 r J 

(33) 

where k = N + 21 — a, with a as the shifting parameter, and U nr (r) is the reduced 
radial wave function. Eq. ( p3| ) is written in units h = c = 1. In these units we have 
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that a = 1/137. Here we borrow the results reported in Ref [25] Introducing the 
scaled mass 

M a = (M 2 + mB) 1/2 (34) 

and proceeding to minimize Eq 

E = V(r ) + (M 2 + Q/(2r ) 2 y/ 2 , (35) 
we obtain that r$ satisfies the equation 

r V(r )(l + Q/(4M a 2 r y/ 2 = Q/AM a . (36) 
The shifting parameter is 

a = 2 - (1 + 2n)w (37) 

where w is given by 

w = ( 3 + r + V"(r )/V'(r ) - 4r V(r ) 2 /Q) ^ , (38) 
and Q can be written as 

Q = k 2 , Q=[r 2 o V'(r )} 2 (2 + 2g), (39) 

with 

g = (l + (2M a /(r V'(r ))) 1 / 2 . (40) 
Equations (^7[) and ( |39"D along with equations (|38| ) and (fH]) give 

2/ + (1 + 2n)w - r 2 V'(r ) (2 + 2#) 1/2 = 0. (41) 

Equation (|41| ) allows computing the value of r . Using this value we compute the 
coefficients 

E x = E V"'(r ) - V(r )V"'(r ) - W(r )V"(r ), (42) 
E 2 = E V""(r ) - V(r )V""(r ) - W \r Q )V m '(r ) - 3y"(r ) 2 , (43) 
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el = (2 - a)/w, e2 = -3(2 - a)/(2w), 



(44) 



£ 3 = -1 + r^x/(3Qu; 3/2 ), e4 = 5/4 + r^ 2 /(12Qw 2 ) (45) 
and finally we have for the energy 

E = E + l/2/£ r 2 ((l - a) (3 - a)/4 + (1 + 2n)e2 + 3(1 + 2n + 2n 2 )e4 (46) 
-l/w(el 2 + 6(1 + 2n)ele3 + (11 + 30n + 30n 2 )e3 2 ) + 5 2 r^/4) 

where the energy level is expressed in units where c = 1, K — 1 . The energy eigenvalues 
can be expressed in atomic units with the help of the relation: 

E n = 2(E-l)/a 2 . (47) 

Now we proceed to compare the results obtained via the mixed variational approach 
as well as the 1/N expansion with those computed numerically. Energy levels are 
computed in Rydberg units as a function of 7' = 2lul/(2lul + 1) and displayed in 
Figures 1 to 4. 

4 Comparison with the 2D nonrelativistic Hydro- 
gen Atom 

In this section we proceed to compare the results obtained for the energy spectrum 
of the 2D relativistic Hydrogen atom with those computed in the nonrelativistic 
limit, when the 2D Schrodinger equation was considered. In order to establish a 
better comparison we use numerical results obtained with the Schwartz interpolation 



method p4|. In Tables 1 and 2 we exhibit different values of the energy for different 



magnetic field strengths. 
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Figure 1: Energy of the 2P~ state as a function of 7'. The solid line is obtained by 
numerical methods; the long-dashed line corresponds to the mixll basis; the short- 
dashed line corresponds to the mix21 basis, (2P~, 3P~ Hydrogen bases and 2P~ 
oscillator wavefunction). The dotted line is obtained by using the mixl2 basis (2P~, 
3P~ oscillator bases and the 2P~ Hydrogen wavefunction). The dash-dotted line is 
obtained with the help of the shifted 1/N method 
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Figure 2: The figure shows the difference between the numeric result for the 2P~ 
energy spectrum and the energy values computed with the help of the mixll varia- 
tional basis (solid line), mix21 variational (light dashed line), mixl2 variational (heavy 
dashed line), and the shifted 1/N method (dotted line) 
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0.2 0.4 0.6 0.8 



Figure 3: Energy of the 3-D - state as a function of 7'. The solid line is obtained by 
numerical methods; the long-dashed line corresponds to the mixll basis; the short- 
dashed line corresponds to the mix21 basis, (3D - , AD~ Hydrogen bases and 3D~ 
oscillator wavefunction). The dotted line is obtained by using the mixl2 basis (3D - , 
AD~ oscillator bases and the 3D - Hydrogen wavefunction). The dash dotted line is 
obtained with the help of the shifted 1/N method 
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Figure 4: The figure shows the difference between the numeric result for the 3D~ 
energy spectrum and the energy values computed with the help of the mixll varia- 
tional basis (solid line), mix21 variational (light dashed line), mixl2 variational (heavy 
dashed line), and the shifted 1/N method (dotted line) 
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7' 


nonrel 


rel 


diflf x 100 





-0.4444444 


-0.44449574 


0.005134 


0.1 


-0.5239504 


-0.52421206 


0.026166 


0.2 


-0.5629132 


-0.5635173 


0.06041 


0.3 


-0.5635941 


-0.56468762 


0.109352 


0.4 


-0.5193465 


-0.52115542 


0.180892 


0.5 


-0.4095808 


-0.41248102 


0.290022 


0.6 


-0.1862031 


-0.19088784 


0.468474 


0.7 


0.2719562 


0.2640041 


0.79521 


0.8 


1.3504374 


1.33516298 


1.527442 


0.9 


5.1012487 


5.0599714 


4.12773 



Table 1 Relativistic energy values for m = — 1, and a comparison with the non relativistic 
energy spectrum. The first column corresponds to the nonrelativistic energy, the second 
column is the relativistic energy E — c 2 , and the third column corresponds to 100 times the 
difference between first and second columns 



7' 


nonrel 


rel 


diff x 100 





-0.16 


-0.1600366 


0.00366 


0.1 


-0.2605089 


-0.2606836 


0.01747 


0.2 


-0.2731927 


-0.27364026 


0.044756 


0.3 


-0.2377636 


-0.2386137 


0.08501 


0.4 


-0.150434 


-0.1518818 


0.14478 


0.5 


0.0113883 


0.00901904 


0.236926 


0.6 


0.3009099 


0.29702154 


0.388836 


0.7 


0.8503051 


0.84361474 


0.669036 


0.8 


2.0734797 


2.06046406 


1.301564 


0.9 


6.1341553 


6.09847572 


3.567958 



Table 2 Relativistic energy values for m = —2, and a comparison with the non relativistic 
energy spectrum. The fist column corresponds to the nonrelativistic energy, the second 
column is the relativistic E — c 2 energy, and the third column corresponds to 100 times the 
difference between first and second columns 

Tables 1 and 2 show that the role played by the relativistic corrections is to shift 
down the energy levels. The relativistic effects becomes noticeable when the magnetic 
field parameter 7' is close to unity. Figure 5. shows the dependence of this difference. 
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3d difference 
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0.2 0.4 0.6 0.8 1 

y/(y+l) 

Figure 5: The dotted line represents the difference between nonrelativistic and rela- 
tivistic 2p energy levels. The solid line represents the difference between nonrelativis- 
tic and relativistic 3d energy levels 
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5 Concluding remarks 



Figures 1 thru 4 show that the variational mixed-basis method gives very good results 
for the energy spectrum of a 2D relativistic Hydrogen atom in a constant magnetic 
field. The energy eigenfunctions obtained with this approach give good energy values 
even for intermediate values of 7' as shown in Figs 1-4 where we compare the varia- 
tional results with those obtained numerically via the Schwartz interpolation method 



23 1 . The shifted 1/N method fails to give reasonably results when we apply its ex- 



tension to the Klein Gordon equation [25[. The advantage of the mixed variational 
method is that we obtain a simple form of the wave function and a reasonable good 
approximation without considering a large term variational basis. The role played by 
relativity consists in shifting down the energy levels as indicated in Figure 5. The re- 
sults presented in this article were obtained considering a two-dimensional hydrogenic 
atom in the presence of a magnetic field perpendicular to the plane of motion. In this 
direction there are some differences between our approach and the method applied in 
Ref. |18 , |i~9f . Here the authors consider a non relativistic quasi two-dimensional sys- 



tem confined by a square- well Vb(z). The Hamiltonian in Ref. |18], [19| contains this 
term and therefore the energy spectrum depends on the dimensions of the confining 
well. Since the variational technique applied in this paper is equivalent to the method 
suggested by Chen et al \TB, HJ, our results reduce, in the non relativisic limit, and 
when the width of the well is negligible, to those reported by Chen et al. Finally we 
mention that s states are not present for the 2D Klein-Gordon Hydrogen atom. The 
absence of s (m = 0) states for the relativistic Klein-Gordon 2D Hydrogen atom can 
be understood if we look at the behavior of Eq. (|TT|) as r approaches to zero. For 
m=0, we have a "falling to center" problem ||26||, and this behavior is unobserved when 
we solve the 2D Dirac equation. A detailed discussion of relativistic effects including 
spin corrections will be presented in a future publication. 
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